Fluctuation theorem in quantum heat conduction 
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We consider steady state heat conduction across a quantum harmonic chain connected to reservoirs 
modelled by infinite collection of oscillators. The heat, Q, flowing across the oscillator in a time 
interval r is a stochastic variable and we study the probability distribution function P(Q). We 
compute the exact generating function of Q at large r and the large deviation function. The 
generating function has a symmetry satisfying the steady state fluctuation theorem without any 
quantum corrections. The distribution P(Q) is nongaussian with clear exponential tails. The effect 
of finite r and nonlinearity is considered in the classical limit through Langevin simulations. We 
also obtain the prediction of quantum heat current fluctuations at low temperatures in clean wires. 
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A lot of interest has been generated recently in fluc- 
tuations in entropy production in nonequilibrium sys- 
tems. Several definitions of entropy production have been 
used and these give some measure of "second law viola- 
tions" . A number of authors have looked, both theoret- 
ically [H, 0, 0, 0] and in experiments [1, 0], at fluctu- 
ations of quantities such as work, power flux, heat ab- 
sorbed, etc. during nonequilibrium processes and these 
have been generically referred to as entropy production. 
The new results, referred to as the fluctuation theorems, 
make general predictions on the probability distribution 
P(S) of the entropy S produced during a nonequilibrium 
process [l|, Q. Specifically these theorems quantify the 
probability of negative entropy producing events which 
become significant if one is looking at small systems or 
at small time intervals. There are two different theorems, 
the transient fluctuation theorem (TFT) and the steady 
state fluctuation theorem (SSFT). The TFT looks at the 
entropy produced in a finite time t in a non-steady state. 
In SSFT, one looks at a nonequilibrium steady state and 
the average entropy production rate over a long time in- 
terval t is examined. The precise statement of SSFT is: 
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In the context of SSFT a quantity of great interest is the 
large deviation function h(a) which specifies the asymp- 
totic form of the distribution function P(S) through the 
relation P(S) ~ e Th( - a) 0,0] • An equivalent statement 
of SSFT can be made in terms of a special symmetry of 
h(a) which is: h{a) — h(—a) = a. Remarkably, this re- 
lation has been shown to lead to linear response results 
such as Onsager reciprocity and the Green-Kubo rela- 
tions @, H, @|. Furthermore it leads to predictions for 
properties in the far from equilibrium regime. 

Heat conduction is a natural example where one 
talks of entropy production. The standard result from 
nonequilibrium thermodynamics is that when an amount 



of heat Q is transferred from a bath at temperature Tjr, to 
a bath at temperature Tr (< Tl) the entropy produced 
S is given by S = (T^ 1 — T£ ) Q. However, in gen- 
eral S is a stochastic variable with a distribution P(S). 
The distribution P(S) for a nonlinear chain connected to 
Nose-Hoover baths was studied numerically in_[10| where 
they verified that it satisfied SSFT. Refn. [llj studied 
heat conduction in a nonlinear chain connected to free 
phonon reservoirs. Based on strong ergodicity properties 
of the model, it was proved that P(S) ~ e rh ^ where 
h{<j) satisfied the SSFT symmetry. In a set-up with di- 
rect tunneling between two finite systems, a transient 
version of the heat exchange fluctuation theorem, valid 
both for classical and quantum systems, was proved in 

While the SSFT clearly presents a powerful theorem 
for nonequilibrium systems, its validity has been estab- 
lished only in specific systems and so far only classically. 
In the transient version, it was proved that quantum 
corrections are necessary for a dragged Brownian par- 
ticle fliSj ]. Thus it is an open question as to whether 
quantum corrections to SSFT exist in quantum heat 
transport and what the characteristics of the heat cur- 
rent distribution are. This letter presents the first ex- 
plicit calculation of h(a) and demonstration of SSFT in 
quantum heat conduction. We study steady state of a 
quantum harmonic chain connected to baths which are 
modelled by infinite oscillator sets. This model is rele- 
vant to recent experiments on mesoscopic quantum heat 
transport fLU HBj ]. where the quantized thermal conduc- 
tance go(T) = irk^T/(6h) was measured 15, 3|. We use 
the method of full-counting statistics [T3| to compute the 
generating function of Q. We then show that the corre- 
sponding large deviation function satisfies the SSFT sym- 
metry condition. For finite r we consider heat transport 
across small chains and study the classical limit through 
Langevin simulations. We also consider the effect of in- 
troducing nonlinearity in the oscillator potential. 
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Model. — Our model consists of a harmonic chain cou- 
pled to two heat baths kept at temperatures Tl and T R 
respectively. For the heat baths we assume the standard 
model of an infinite collection of oscillators. The full 
Hamiltonian is given by 
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where {m n ,x n ,Pn,kn,k} refer to the system degrees of 
freedom, , m^, uji\ refers to the left reservoir while 

{x r ,p r ,m r ,u) r } refers to the right reservoir. The cou- 
pling constants between the system and the bath oscil- 
lators {A^,A r } is switched on at time t = — oo. The 
initial density matrix is assumed to be of the prod- 
uct form p{— oo) = ps ® Pl ® PR,, where S,L,R re- 
fer respectively to the system and left and right reser- 
voirs. The left and right density matrices are equilib- 
rium distributions corresponding to the respective tem- 
peratures: p a — e~ /3o,Wa /Tr[e~ /3o,Wo! ] for a = L,R and 
P a = l/(fc B T Q ). 



It can be shown [l8| that eliminating the bath degrees 
of freedom leads to an effective quantum Langevin equa- 
tion for the system. The effect of the baths is to produce 
noise, given by r]L jR (t), and dissipative effects controlled 
by memory kernels jL,Fi(t). The properties of the noise 
and dissipation are completely determined by the initial 
condition of the baths at t — — oo. We now make a few 
definitions. Let rj LtR (u}) = dtj L>R (t)e wt , Jl,r(u)) = 
J °° dt^ L , R (t)e luJt and let Y7 l r {lj) = -iuj^ L , R (oj), which, 
as we will see later, gives the self energy correction 
coming from the baths to the Green's function of the 
harmonic chain. We also define the spectral function 

A 2 

Jl(u) — f ^2 e S(ui — uj() for the left reservoir and a 
similar function J r (lo) for the right reservoir. Then the 
dissipation kernels and noise correlations are given by: 

2 r°° j q h 

7 a it) = — I dui cos ut 

faaMfafw')) = 47rM(a; + w')r a (w)(l + / a (w))(3) 



for a — L,R and where T a (uj) = — Im{S^(ijj)} = 
J a (u)Q(u) - J a (-w)Q(-(jj) and f a (w) = l/( e ^ huJ - 1). 
All higher noise correlations can be obtained from the 
two-point correlator. Using the quantum Langevin ap- 
proach it is straightforward to derive the Landauer type 



result for average heat current (X) [LSI ]: 
(I) = -L / dujhujT(uj)[ f L (u) - Jr{u)) 



(4) 



T(w) = 4T l (uj)T r (ll>)\G t 1n (uj)\ 2 , 

G» = [Muj 2 — K - E£(w) - VUoj)}- 1 , (5) 

where M, K are the mass and the force constant matrix 
and H r L r (lo) are self-energy correction matrices with el- 
ements [££(w)] m ,„ = E£(w)£ m ,n<W and [XJ^(w)] m , n = 
Ti r R {uj)5 m ,n5m,N ■ Note that T(w) is the transmission co- 
efficient for phonons while G r (oj) is the phonon Green's 
function for the chain. 

Statistics of phononic heat transfer. — The heat trans- 
fer in time r is given by (T)t. Here we are interested 
in the statistics of heat transfer in the nonequilibrium 
steady state and so need to calculate higher moments 
of the heat transfer. The quantum Langevin approach 
can in principle be used to compute correlation functions 
at any order but this becomes increasingly cumbersome. 
Instead we use the Keldysh approach which, as we will 
show, gives the generating function of the heat transfer. 

Several definitions of X are possible depending on 
where we evaluate the current. Here we consider the 
current from the left reservoir into the system ( obtained 
by taking a time-derivative of energy in the left reservoir 
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We also define the average heat transfer operator Q = 
J-t/2 dtT(t) • Using the Keldysh approach let us compute 
the following quantity: 
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where (...) denotes an average over the initial state, T 

and T denote forward and reverse time ordering, and 
the counting field <p(t) is defined as <p(t) — —h^/2 for — 
t/2 < t < t/2 and zero elsewhere. It can be shown that 
lnZ(£) is the cumulant generating function for the heat 
operator: 



lnZ(0 = 



71=1 
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(7) 



where {Q n ) c is the n order cumulant at large r. 
Hence the probability distribution of measuring a heat 
transfer Q is obtained by taking the Fourier transform 
P(Q) = ± d£ Z(£) e~^ Q . For large r one obtains 

Z(0 ~ e rG ^) and P{Q) ~ e T ~ h ^ with h(q) = G(C)-i£*q 
and where £* is the solution of the saddle-point equation 
dQ(£*)/d£* — iq = 0. We evaluate <?(£) using standard 
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path integral and Green's function techniques along the 
Keldysh contour. The final result is the following form: 

-1 f°° ( 

S(0 = ^ J ^ln{l+TH[/ fl (- W )/ L H(e^-l) 

+ f R (u;)f L (-u)(e-^-l)}} . (8) 

Phonons convey energy in units fko and this appears in 
the exponential form with the factor £. It is easily verified 
that Eq. ijHJ reproduces the correct first moment of X 
given in Eq. The second moment is given by 

^ = ^ /°° (M 2 {^M [/i(«)-Wf 

- r(u)[/ i (w)/jiH+/ i (-w)/ J e(u)]}. (9) 

We have verified this also with the Langevin approach. 
This bosonic fluctuation is similar to the optical one • 
Symmetry. — We note the following symmetry of Q: 

G(0=g(-Z + iA) , (10) 

where A = Pr — (3l- Using the identification a = Aq and 
the relation between h(a) = h(q) and £/(£) immediately 
leads to the SSFT relation Eq. ([T]). Thus we conclude 
that quantum heat transports satisfy the SSFT without 
any quantum corrections. 

The symmetry (|10[) contains information regarding 
transport coefficients J§|. For fixed (3l + (3r let us 
make the expansion (<2") c /t = J2 m L n , m A m /m\. The 
nonlinear response coefficients L n>m are then given by 
L n>m = d n+m g{£)/d(i£) n dA m \ i=A=0 . This coefficient 
represents a nonlinear response of the general cummu- 
lants of current to the thermodynamic force (/3r — Pl)- 
The symmetry (|TU|) gives the general relations between 
the coefficients: 

L n>m = J2 (™) (-l) {n+k) L n+k , m -k, (11) 

with Lo, m = 0. For example we get L2.0 — 2Li i and 
L4 = 2L31 = 6L2.2 — 4:Li,3- The first relation is simply 
the Green-Kubo formula relating the linear current re- 
sponse to equilibrium fluctuations while the second leads 
to relations between nonlinear response coefficients. 

Typical distributions. — We present some results on the 
form of the distribution P(Q) ~ e h ^ for a small chain 
(N = 2) connected to ohmic reservoirs {^l,r{u) = 7)- In 
Fig.[T]we plot h(q), which is numerically obtained, for dif- 
ferent temperatures Tl with fixed temperature difference 
Tr — Tl- In all temperature regimes, h(q) shows a clear 
linear dependence at large q, and those are well fitted by 
0Rq and — for g < and q > respectively. This 
exponential tail is one of the characteristics in P(Q). 

We now study the effects of a finite r and nonlinear 
potential using the classical system. We evaluate P(Q) 
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FIG. 1: Plot of h(q) for various temperature regime with Tr — 
Tl fixed to 1.0 [hk^ y/k/m\. The parameters:^ = 2, mi = 
mj = m, ki = k-i = k, and 7 = 1.0 [vmfe]- 




FIG. 2: Plot of Ih[P(Q)]/t from Langevin simulations for the 
same system as in Fig. [T]for various r. Here P(Q) is normal- 
ized so that its maximum is one. The large deviation function 
h(q) is also shown by black solid line. The parameters are 
N = 2,m = k = 1, k B T L = 10, and k B T R = 20. The non- 
linear case has a = 2. The inset shows Ih[P(Q)/P(—Q)]/t, 
and SSFT line q(/3n — (3l) (black solid line). Heat transfer is 
measured at the contact to the left reservoir. 

from direct simulations of the classical Langevin equa- 
tions with white noise. In Fig. [5] we compare the simula- 
tion results for different values of r with the asymptotic 
function h(q) (obtained for H — ► 0). It is clear from Fig. [2] 
that convergence to the asymptotic distribution function 
takes place on a rather large time scale. The nonlinear 
case is also plotted for the same system with an onsite 
potential V(x n ) — In the inset, the function 

ln[P(Q)/P(— Q)]/t is plotted for three cases. The distri- 
bution for the nonlinear case deviates from the harmonic 
cases, and both the average heat current and its fluc- 
tuations are suppressed. However, as the inset shows, 
SSFT is satisfied in the nonlinear case, which indicates 
the symmetry (fTU)) and the relation (jTTj) hold too. 

Heat current fluctuations in a pure wire. — Using 
Eq.®, we can derive the heat current fluctuations for a 
homogeneous wire connected to reservoirs through non- 
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reflecting contacts, a case for which the quantized ther- 
mal conductance has been measured [l5|. Consider a 
pure wire with all masses and spring constants equal. If 
we consider that the heat reservoirs themselves are semi- 
infinite wires (i.e., the Rubin model of a heat bath) then 
it is easy to show that the contacts are perfect and we get 
T(u)) = 1 for all lu within the allowed bandwidth. At low 
temperatures and for small AT = Tl — Tr, Eq. (J4j) leads 
to the quantized heat conductance go(T) = (2) /AT = 
Trk%T/(6H). From Eq. ((9]) we now also get the thermal 
noise power at zero frequency Sq = (Q 2 )/t: 
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Sq = k B Tlg (T L ) + k B T^g (T R ). 



(12) 



This is valid for Tl,_r in the temperature regime where 
go(T) can be measured [15j . The noise power is also inde- 
pendent of details of system. Independent contributions 
from Tl and Tr are obtained since there are no scattering 
process between phonons. Eq. (jSJ) with T(u) = 1 gives us 
the generating function valid in the same regime. While 
our results have been derived for a one-dimensional wire 
with scalar displacement variables they are easy to gener- 
alize. Similar results can be obtained for realistic models 
[20| of nanowires and nanotubes. 

Summary. — Unlike equilibrium physics there are few 
general principles to describe nonequilibrium phenom- 
ena. The exceptions to this are the Onsager reciprocity 
and the Green-Kubo relations which are valid in the 
close-to-equilibrium linear response regime. In view of 
this the nonequilibrium fluctuation theorems are quite 
remarkable in that they seem to be exact relations valid 
arbitrarily far from equilibrium and from which one can 
recover standard linear response theory. However the 
full range of validity and applicability of these theorems 
is still not known. In this paper we have derived the ex- 
plicit distribution for fluctuations in phononic heat trans- 
fer across a quantum harmonic chain and have obtained 
the first proof of SSFT in quantum heat conduction. We 
find that there are no quantum corrections. We note that 
fluctuations in charge current in mesoscopic systems have 
been already studied both theoretically [l?], O, 21, [22| 
and experimentally [23] |. The present study provides a 
theoretical basis to study flucutuation of heat transfer. 
The measurement of fluctuations of phononic heat trans- 
fer in experiments is an important challenging problem. 
A modification of the set up used in [l5| should be able 
to make a measurement of fluctuations in heat trans- 
fer. One possibility is to use heaters with some feed-back 
mechanism so as to maintain the two reservoirs at fixed 
temperatures. The fluctuations in the power from the 
heater would be related to the fluctuations in the heat 
transfer through the wire. 
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